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Review  of  Fourier,  Lapiecs  and  Z-transforms  provides  introduction  to  discrete  time  systems. 
The  KC  low-pass  filter  is  modeled  in  terms  of  a  feedback  structure  and  compared  with  its  analogy 
to  recursive  low-pass  digital  filters.  These  topics  introduce  the  preferred  transversal  digital 
filter  for  low-pass  and  band-pass  applications.  An  efficient  method  of  synthesizing  transversal 
filters  is  demonstrated  from  the  derivation  of  Fourier  coefficients  and  time  delays.  A  method 
of  weighting  the  truncated  Fourier  terms  by  Hamming  coefficients  to  suppress  Gibb’s  oscillation 
'  is  also  shown.  Design  examples  demonstrating  the  magnitude  and  phase  response  of  low-pass 
and  band-pass  transversal  filters  are  presented. 
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INTRODUCTION 


The  rise  of  digital  filtering  has  come  about  as  a  result  of  two  factors  which 
are  not  totally  independent.  The  first  is  traceable  to  the  fact  that  more  process 
cc  ,rol  and  signal  processing  operations  are  now  being  executed  digitally  which 
were  previously  accomplished  in  the  analog  domain.  Once  signals  are  in  a  digital 
format,  it  is  obviously  more  efficient  to  complete  all  necessary  operations  digitally. 
The  second  factor  bringing  about  a  more  widespread  application  of  digital  filters 
is  the  continuing  decrecse  in  cost  per  digital  hardware  function.  This  is  due  to 
<•’  'vances  in  the  state-of-the-art  and  economic  factors  of  high  volume  production. 
Digital  filters,  therefore,  may  be  realized  with  either  software  as  a  part  of  a  larger 

f 

executive  program,  or  in  a  hardwired  hardware  form,  with  perhaps  A/D  and  D/A 
input  and  output  converters,  to  accomplish  a  specific  filtering  task. 

The  decision  of  when  to  use  digital  filters  instead  of,  for  example,  active- 
type  analog  filters  is  not  dictated  by  hard-and-fast  guidelines.  This  decision  is 
strictly  dependent  upon  the  application  at  hand.  The  principal  features  of  digital 
filters  are  high  stability  and  the  ease  with  which  filter  parameters  can  be  changed. 
The  former  is  useful  for  separating  signals  very  close  in  frequency  and  where  the 
filter  must  be  located  in  other  than  ideal  environments.  The  latter  feature  is  essential 
for  adaptive  filtering  where  a  filter  must  change  its  characteristics  in  real-time  to 
a  changing  input  signal.  The  following  development  presents  the  essentials  of 
digital  filtering  from  a  practical  viewpoint  which  the  practicing  engineer  can 
immediately  put  to  use  and  includes  both  low-pass  and  band-pass  design  examples 
with  their  output  response  plots. 


FOURIER  TRANSFORMS 


The  Fourier  transform  pair  is  useful  for  transforming  functions  between  the 
time  and  frequency  domains 

g(“)  -  /  Kt).'1"'* 

-CD 

f(»)  *  -jjr  /  fl(«)  •  d« 

-CO 

Consider  a  single  pulse  of  amplitude  A  and  width  t.  Its  transform  from  the  time 
to  frequency  domains  provides  the  spectral  density  of  the  pulse.  It  can  be  seen 
that  a  pulse  of  width  t  requires  a  bandwidth  of  2/t  Hertz  for  transmission  or  re¬ 
production  with  fidelity  since  the  majority  of  the  energy  contained  within  the  pulse 
is  seen  to  lie  between  +  1/t  Hertz  of  the  spectral  density  plot. 


T  T 
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LAPLACE  TRANSFORMS 


The  Laplace  transform  pair  Is  a  one-sided  Fourier  transform  pair  useful  for 
transforming  discontinuous  functions  between  the  time  and  frequency  domains  which 
the  Fourier  transforms  cannot  handle.  For  example/  consider  a  Fourier  transform 
of  a  unit  step  function 
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This  result  is  undefined  and  therefore  meaningless.  However,  convergence  of  this 
integral  can  be  assured  by  adding  a  constant  a  to  jw  providing  s  =  a  +  ju  ,  which 
is  a  complex  variable.  Considering  now  the  Laplace  transform  of  a  unit  step 
function  we  obtain  a  useful  solution 

g(a  +  ju)  =  /  1  *  e  dt 
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Z  -  TRANSFORMS 


Laplace  transforms  are  principally  useful  for  continuous  functions.  If  the 
Laplace  transform  of  a  time-sampled  continuous  function  is  considered 


F(s)=/  f(t)  e  dt 
o 


n=0 


The  final  expression  evaluates  each  time  t  =  nT.  For  convenience  we  defi.ie  a  unit 
time  delay  operator  z  ^  =  e  which  provides  the  Z-transform  pair 

F(z)  =  I  f(nT)Z“n 

n=0 

f(nT,=-l_  /  F(z)Zn-1dz 
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RC  LOW  PASS  FILTER 


Both  an  analog  and  digital  realization  of  a  single-section  RC  filter  will  be 


considered  as  an  illustration  of  an  intuitive  approach  to  digital  filter  design. 


TM  -  -  1 

T  1  +  RCS 


Recalling  that  the  transfer  function  for  a  feedback  structure  is  given  by  T  =  G/O+GH), 
we  represent  the  analog  RC  filter  in  like  manner 


Consider  now  filter  excitation  by  a  unit  step  function  beginning  at  t  =  0.  The  result 
ing  RC  filter  output  response  is  obtained  with  the  aid  of  Laplace  transforms  of  this 
input  step  and  the  filter  transfer  function 


X(s)  =  /  1  •  e  dt 
o 


Expansion  of  Y(s)  by  partial  fractions  and  obtaining  the  inverse  Laplace  transform 
from  tables  yields  the  RC  filter  output  tir  j  response  in  terms  of  an  exponential 
function 


Y(t)  =<£.  1Y(s)1 

=  l-e"^C 
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RECURSIVE  DIGITAL  FILTERS 


The  foregoing  single-section  RC  filter  may  also  be  realized  in  digital  form 
by  obtaining  the  Z-transform  of  the  input  step  function  and  the  just  presented  output 
time  response  Y(t).  The  transfer  function  formed  with  these  quantities  is  also  realiz¬ 
able  with  a  recursive  feedback  structure  as  shown  for  the  analog  case. 


X(z)  =  T  X(nT)Z 
n=0 


The  foregoing  introduction  Is  meant  to  bridge  the  gap  between  continuous- 
function  analog  filters  and  sampled-function  digital  filters.  The  particular  example 
chosen  is  intended  as  a  simple  exercise  rather  than  a  practical  example  of  digital 
filter  design.  Recursive  filter  designs  are  feasible  only  under  rather  severe  realizability 
constraints.  The  details  of  these  constraints  can  only  be  briefly  addressed  here. 
However,  since  they  are  feedback  -'♦ruJures  recursive  filters  are  subject  to  all  of 
the  stability  considerations  of  analog  feedback  systems  plus  additional  ones  which 
can  defy  practical  solution.  For  example,  quantization  noise  encountered  in  the 
sampling  process  when  fedback  and  multiplied  by  the  various  filter  coefficients  can 
produce  non-linear  behavior  and  filter  limit  cycling.  In  addition,  passband  phase 
shift  in  recursive  filters  is  generally  difficult  to  linearize  and/or  compensate  satis¬ 
factorily.  A  more  practical  and  useful  approach  to  digital  filters  is  possible  with 
the  transversal  structure.  For  the  above  stated  reasons  the  following  development 
is  devoted  to  a  detailed  presentation  of  transversal  digital  filter  design. 
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TRANSVERSAL  DIGITAL  FILTERS 


The  transversal  filter,  also  known  as  a  finite-memory  or  nonrecurstve  filter, 

Is  not  suitable  for  Integration  applications  due  to  this  availability  of  only  a  finite 
memory.  However,  it  Is  useful  for  low-pass  filtering,  interpolation  and  extrapolation 
(delay  and  prediction),  differentiation  and  combinations  of  the  above.  In  addition, 
and  to  be  introduced  subsequently,  with  the  low-pass  to  band-pass  transformation  the 
transversal  filtermechanization  provides  a  band-pass  filter  with  all  of  the  advantages 
represented  by  transversal  filters. 

Digital  filtering  of  continuous  signals  involves  the  signal  processing  operations 
of  sampling  and  quantization,  which  is  addressed  elsewhere  in  this  treatise,  as  well 
as  the  computational  operations  of  storage,  multiplication  and  addition.  Digital 
filters  can  be  ma-Je  to  approximate  their  RLC  counterparts  insofar  as  conforming  to 
the  Butterworth,  Bessel,  or  Chebyshev  characteristics,  for  example,  if  this  is  of  in¬ 
terest.  However,  a  more  direct  approach  is  elected  in  this  presentation  whereby 
we  will  be  concerned  with  considerations  of  passband,  transition  regions,  stopband, 
ripple  and  phase.  This  direct  approach  will  provide  more  insight  into  the  design  of 
digital  filters  without  manipulating  esoteric  transfer  functions  in  order  to  achieve 
particular  response  characteristics. 

Low-pass  transversal  digital  filter  structures  are  useful  for  estimating,  at  discrete 
time  instants,  the  value  of  a  random  signal  perturbed  by  additive  noise.  Transversal 
filters  are  always  stable  since  no  feedback  paths  are  involved  and  designs  with  exactly 
linear  phase  are  easily  realized.  In  addition,  due  to  this  inherent  stability  coefficient 
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value  approximations  may  be  relaxed  which  permits  shorter  coefficient  multiplier  word- 
lengths.  This  plus  the  fad  that  only  a  finite  number  of  states  are  generally  required 
for  a  transversal  filter  realization  provides  efficiency  ?n  its  hardware  mechanization. 
The  spectral  response  in  the  Z-domain  of  a  transversal  filter  is  expressed  by  Equation 
(1)  and  its  block  diagram  by  Figure  1 . 

H(z)  =  I  a  Z"n  (1) 

n=0 


X(nT) 


Figure  1.  Transversal  Filter  Mechanization 


LOW  PASS  TRANSVERSAL  FILTERS 


Consider  now  the  design  of  a  low  pass  filter  whose  frequency  response  is  in 

general  represented  by  Figure  2.  In  radians  per  second,  u  is  the  filter  cutoff 

c 

frequency,  the  frequency  where  high  attenuation  exists,  and  <*>s  the  filter 
sampling,  computation  and  Z  ^  delay  rate  which  in  practice  is  selected  sufficiently 
high  to  satisfy  the  Nyquist  Sampling  Theorem  and  to  prevent  overlap  of  the  periodic 
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spectrum  foldover.  This  will  be  addressed  in  sufficient  detail  at  the  appropriate 
point  later. 


Figure  2 .  Frequency  Domain  Filter  Response 

The  designer  initially  determines  the  allowable  filter  transition  region  between 
the  passband  and  the  stopband  which  is  given  by  Equation  (2) 


Since  there  are  no  signal  components  passed  by  the  filter  at  frequencies  above  u 

stop 

then  the  filter  output  must  drop  from  its  passband  value  at  u  to  the  desired  attenuation 

c 

at  “stop*  F'lter  ro,,off  terms  of  DB/decade  or  DB/octave  is  therefore  easily  specified. 
Note  that  no  comment  was  made  about  how  many  DB  down  the  filter  respcr.se  will  be 
at  the  cutoff  frequency.  This  is  because  it  is  not  necessary  to  impose  the  parametric 
characteristics  of  analog  Rl£  filters  on  digital  filter  designs.  As  we  will  see  shortly, 
response  factors  are  determined  by  parameters  alien  to  analog  filter  design.  Before 
discussing  this,  however,  the  Fourier-series  method  for  obtaining  the  filter  multiplier 
coefficients,  a^  of  Figure  1,  will  be  presented. 


In  practice  ,  the  series  is  truncated  symmetrically  for  economy  of  implementation. 

This  yields  a  total  of  2J+1  filter  states,  where  m=2J.  The  Z-domain  expression  for 
the  transversal  filter  in  terms  of  Fourier  coefficients  is  given  by  Equations  (10)  and  (11) 
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(10) 


J  J 

H(z)  =  a  i-m  (a  -ft)Z'n+I  (a  +ft)Zn) 


o  ,  n  '  n 

n=l 


=  1  a  r 

n=0  .. 


The  filter  coefficient  multipliers  are  specifically  given  by 


•2  (  Oj  +  P  j) 

!  ^J-l  +  PJ-1} 

=  •2  (a,  +  P,) 

•  i 

-  a 

o 

=  2  (a^  “  (Jj) 


Vi  =*  <0J-i  ?j-i* 

°m 

If  the  Fourier  series  is  truncated  using  a  rectangular  window  as  indicated  by 
Equations  (10)  and  (11)  and  the  resu'  ing  filter  multiplier  coefficients,  the  filter 
frequency  response  will  exhibit  Gibb1  s  oscillation  due  to  the  discontinuity  in  H(u). 
Gibb's  oscillation  can  be  effectively  reduced  by  modifying  the  Fourier  coefficients 
by  means  of  a  weighting  function.  The  Hamming  window  is  particularly  useful  for 
controlling  the  ripple  in  the  frequency  response  of  finite-length  sequences.  A  function 
describing  this  window  is  presented  by  Equation  (12) 

0.54  +  0.46  cos  ( ^-5-)  n  <  J 

\A/  _  i  J 


n  >  J 
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Use  of  the  Hamming  weighting  function  provides  not  only  control  of  the  ripple /  but 
specifies  the  optimum  truncated  filter  sequence  length  J  in  terms  of  the  desired  filter 
transition  region  A  u.  Sampling  time  T  is  defined  in  terms  of  the  sampling  frequency 
where  T  =  2ir/us. 


The  Hamming  weighted  filter  coefficient  multipliers  are  then  given  by 

ao  =  2VV(J)(aJ  +  P|) 

a1  =  *  W(J-1)  (aJ-1  + 

•  • 

•  • 

•  • 

aJ~l  =  -  W(l)  (°1  +  M 

a  ,  =  a 
J  o 

aJ+i  =  2  W(i)  (al  " 

•  • 

•  • 

•  • 

=  =  W(J)  (  °  J  - 

The  Fourier  coefficients  are  obtained  by  noting  that  the  filter  response  in  the  frequency 
domain  is  periodic  at  a  foidover  repetition  frequency  defined  by  the  sampling  frequency 
us>  From  Figure  3  it  can  be  seen  that  in  order  to  prevent  overlap  of  the  periodic 
spectrums,  u$  must  be  sufficiently  high  to  assure  adequate  spacing.  A  rule-of-thumb 
is  that  uc  is  taken  at  two  to  ten  times  the  highest  frequency  of  interest.  If  overlap  is 
allowed  to  occur,  the  phenomenon  of  aliasing  will  take  place  producing  errors  in  the 
reconstruction  of  a  continuous  signal  at  the  output  from  its  sampled  representation. 
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Figure  3.  Periodic  Spectrum  Foldover 

Due  to  the  even  symmetry  imposed  by  the  spectral  response  of  H(u)  about  zero  for  the 

low-pass  filter,  all  of  the  B  terms  will  be  zero.  This  leaves  the  DC  and  even  a 

rn  n 

terms  to  be  evaluated.  Considering  u/2  =  rr/T  we  then  write 


+u 


.  /  n  2ir  u 


pn  =  —  S  c  A  •  Sin  ( ' 
n  u  u, 

S  ~u  S 
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)  d  u  -  0 


04) 


2u 


°o  u 


_  2  a  c  .  » n  2it  U'  \  , 

a  _  —  f  A*cos(  —  )  d  u 

n  u$  u$ 

uc 


(15) 

(16) 


Evaluating  these  expressions  in  the  frequency  domain  with  a  normal lized  amplitude 
A  =  1  yields 


=  —  S  C  1  *  cos  (n  u  T)  d  io 


n  to 


s  -u 


u 


■?  t 


c  1 


nT 


n  T  •  cos  (n«T)  d  w 


=  —  sin  (n  w  T) 

nir  ' 


(17) 


LOW  PASS  DESIGN  EXAMPLE 


Consider  a  low-pass  filter  cutoff  frequency  of  1Hz,  or  u  -  2tr.  Such  a  filter 

could  be  applicable  to  thermocouple  inputs  at  a  data  logging  station.  A  sample 

frequency  must  then  be  chosen  to  be  greater  than  twice  the  highest  frequency  of 

interest  u  .  .  For  a  selected  transition  region  of  Au  =  2ir  the  filter  has  aru,  - 

stop  stop 

4rr.  A  value  for  =  16ir  is  then  selected  as  two  times  the  minimum  practical  sampling 
frequency  for  this  fi  Iter.  Then,  T  =  2ir/us  =  1/8  second  and  the  number  of  required 
filter  states  to  realize  this  filter  is  obtained  from  Equations  (18)  and  (19).  The  Fourier 
DC  and  even  coefficients  are  found,  respectively,  by  Equations  (20)  and  (21) 


J  =  =  8 

T  A  u 

(18) 

m  =  2J  =  16 

(19) 

2uc 

a  ~  C  =  1/4 
o  u 

(20) 

s 

_  2  .  ,  nrr . 

a  -  —  •  sin  ( -T-) 
n  nir  4 

(21) 

The  Fourier  and  Hamming  coefficients  for  this  filter  are  then  obtained  from  Equations 
(21)  and  (16)  as 


FOURIER 

HAMMING 

a  i  =  +  0.462 

WH1  «  +1 .44 

a2  =  +0.318 

=  +  1.16 

Qg  =  +  0.150 

WH3=  +  0.76 

o 

o 

ii 

-^r 

a 

=  +  0.25 
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The  17  filter  multiplier  coefficients  weighted  by  the  rectangular  (simple  truncation) 
and  Hamming  windows  are 
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The  magnitude  and  phase  response  of  the  low-pass  digital  filters  implemented 
with  the  preceding  Hamming  and  rectangularly  weighted  multiplier  coefficients  were 
obtained  from  the  following  frequency-domain  equations 


HU  =  ?  a  e-i"“T 

n=0  n 


=  I  a  (cos(nuT)  -  i  sin  (n«T)l 
n=0  n 


O 


s?n(uT) 
cos(u  T) 


(22) 

(23) 


A  computer  generated  plotting  routine  was  executed  with  these  eauations  for  T  =  1/8 
with  u  ranging  from  0  to  4rr  and  m  =  16.  Figures  4  through  7  present  these  results. 

Consider  the  Hamming  weighted  filter  of  Figure  4.  The  attenuation  achieved 
between  the  passband  magnitude  of  1 .47  and  the  stopband  ripple  peak  of  0.15  is  a 
ratio  of  10,  or  20  DB.  The  transition  region,  however,  extends  over  1 .25  Her*z 


instead  of  the  predicted  1  Hertz.  In  addition,  the  cutoff  frequency  of  1  Hertz  is 
centered  about  halfway  along  the  transition  region  instead  of  at  the  beginning  as 
expected.  These  anomalies  are  primarily  due  to  an  insufficient  sampling  rate. 

Referring  to  Equation  (18),  it  can  be  seen  that  for  the  specified  transition  region 
of  1  Hertz  that  a  higher  sampling  rate,  implicit  in  a  smaller  T,  will  provide  additional 
Fourier  coefficients  and  hence  filter  states  to  better  approximate  the  desired  response. 
The  phase  response  of  this  filter,  presented  by  Figure  5,  is  linear  and  well  behaved 
even  at  this  minimum  sampling  rate.  The  magnitude  and  phase  response  for  the  digital 
filter  formulated  from  rectangularly  truncated  Fourier  coefficients  is  shown  by  Figures 
6  and  7.  The  behavior  of  this  filter  departs  from  the  Hamming-weighted  filter  primarily 
in  terms  of  the  undesired  Gibb' s  oscillation  it  exhibits  in  the  passband.  Illustration 
of  this  oscillation,  of  course,  is  the  purpose  for  providing  this  comparison. 


Reproduced  from  Figure  4.  Hamming  Weighted  Filter  Magnitude  Response 

best  available  copy. 


Figure  5.  Hamming W< 


Figure  6.  Rectangularly  Weighted  Filter  Magnitude  Response 


Figure  7.  Rectangularly  Weighted  Filter  Phase  Response 


BAND  PASS  TRANSVERSAL  FILTERS 


A  low-pass  to  band-pass  digital  ft  Iter  transformation  may  be  effected  by  a 

technique  which  parallels  that  of  a  heterodyning  process  in  analog  signal  theory. 

Performing  this  transformation  on  transversal  filters  allows  the  stability  and  well- 

behaved  characteristics  of  this  mechanization  to  be  realized  for  the  band-pass  case. 

This  straightforward  transformation  is  achieved  in  the  Z-domain  by  the  following 

algorithm  which  mathematically  converts  the  low-pass  filter  multiplier  coefficients, 

a  ,  to  band-pass  multiplier  coefficients,  A  . 
n  n 
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=  Z 
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a  *  cob  [u  T(n-2r)  | 
n  -  r  o 


C27) 


where  u  is  the  filter  center  frequency.  Derivation  of  the  band-pass  coefficients 
o 

from  Equation  (27)  requires  close  attention  to  the  subscripts  on  the  coefficient  an_r* 

Note  that  the  summation  of  terms  on  r  for  conversion  of  the  low-pass  filter  coefficients, 

a  ,  to  the  band-pass  coefficients,  A  ,  is  continued  only  until  n  =  r.  Negative  sub- 
n  n 

script  values  are  thereby  precluded. 
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BAND  PASS  DESIGN  EXA^'lE 


Consider  a  band-pass  filter  based  on  the  preceding  low-pass  transversal  design 
with  the  purpose  of  providing  DC  attenuation  while  passing  slowly  varying  AC  components. 
A  center  frequency  of  2  Hertz  is  chosen  in  the  example  at  hand.  Equation  (28)  provides 
coefficient  conversion  of  the  previous  17  low-pass  Hamming-weighted  coefficients 
with  the  same  sampling  time  of  T  =  1/8  second. 

A  =  a  •  cos  [ir/2  (n-2r)  )  (28) 

n  r=0  n"r 
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Figures  8  and  9  present  the  magnitude  and  phase  response  also  generated  by 
computer  methods  utilizing  Equations  (22)  and  (23)  as  before.  Both  passband  and  stop- 
band  ripple  is  in  evidence  even  though  conversion  was  effected  from  Hamming-weighted 
coefficients.  This  anomaly  is  characteristic  of  this  method  of  generating  a  band-pass 
characteristic.  Ripple  reduction  to  a  specific  acceptable  level  is  achieved  as  a  trade¬ 
off  with  filter  sampling  rate  in  the  same  manner  as  for  the  low-pass  design.  Again, 
from  Equation  (18)  it  is  apparent  that  the  faster  the  scmpling  rate  the  greater  the  number 
of  filter  states,  hence  coefficients  and  delays,  and  thereby  the  more  closely  the  filter 
realization  will  approximate  the  desired  characteristic.  In  the  example  at  hand  12  DB 
attenuation  is  achieved  at  DC.  However,  it  is  apparent  that  a  judicious  selection  of 


ed  F 


gure  9.  Hamming  Weighted  Filter  Phase  Response 


filter  center  frequency/  such  as  1 .6  Hertz,  will  place  the  stopband  null  of  Figure  8 
at  DC  for  24  DB  rejection.  The  phase  response  in  the  filter  passband  i  found  to  be 
reasonably  linear  as  shown  by  Figure  9.  The  phase  warp  shown  occurs  innocuously 
at  the  nulls  of  the  stopband . 
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